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In this paper we study the dynamics of self gravitating spherically symmetric thin 

shells of counter rotating particles. We consider all possible velocity distributions for 

the particles, and show that the equations of motion by themselves do not constrain 

j^j ' this distribution. We therefore consider the dynamical stability of the resulting 

^> ■ configurations under several possible processes. This include the stability of static 
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configurations as a whole, where we find a lower bound for the compactness of the 
shells. We analyze also the stability of the single particle orbits, and find conditions 
for "single particle evaporation" . Finally, in the case of a shell with particles whose 
angular momenta are restricted to two values, we consider the conditions for stability 
£NJ ' under splitting into two separate shells. This analysis leads to the conclusion that 

under certain conditions, that are given explicitly, an evolving shell may split into one 
or more separate shells. We provide explicit examples to illustrate this phenomenon. 
We also include a derivation of the thick to thin shell limit for an Einstein shell 
(3JT)! that shows that the limiting distribution of angular momenta is unique, covering 

continuously a finite range of values. 
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Q , Constructing explicit solutions of Einstein's equations with matter that satisfies a phys- 

ically acceptable equation of state, or, even acceptable energy conditions, has proved to 
be in general a quite difficult task. The problem may be simplified by the introduction 
of symmetries, for instance assuming spherical symmetry, and the imposition of extreme 
equations of state, such as, e.g., constant energy density. A further important geometri- 
cal simplification concerns the assumption that the support of the region where the matter 
density is non vanishing is an embedded hypersurface in the space time manifold, that is, 
that the matter is restricted to a shell. In some applications, such as the recent brane 
models, this is an essential assumption of the theory. In other cases, such as those where 
astrophysical applications are in mind, they should be considered as simplifying limits, and 
therefore, it is quite important to establish what are the configurations for which the shell 
is an idealized limit, and whether the resulting configurations are stable. As regards the 
simplifications on the type of matter, a particularly interesting class, because of its relevance 
in astrophysics, is the assumption that it is made out of collisionless particles, where each 
particle follows a geodesic of the mean field generated, at least in part, by the rest of the 
particles, and where the evolution of the gravitational and matter fields is governed by the 
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Einstein - Vlasov equations [l[. The problem here is again that the general equations turn 
out to be complicated, and one must resort to further simplifying assumptions in order to 
get definite answers. This has led to the consideration of spherically symmetric models, 
where the particle trajectories are also restricted. Two of these restrictions have proved to 
be particularly useful. The first is the restriction to radial motions in the dust models. The 
other corresponds to the assumption that at each point there is a radially comoving frame, 
where the radial velocities of all particles vanish. In the case of spherical symmetry, this 
implies that for each particle moving in a certain (tangential) direction in the comoving 
frame, there must be another particle of the same mass, moving in the opposite (tangential) 
direction. This configuration was first considered by Einstein, and is known as one with 
counter rotating particles. If the particles are further restricted, at least in some spacelike 
hypersurface, to a region R in < r < R out , where r is the surface area radius in the spherically 
symmetric space time, we denote the system as an Einstein thick shell [2L and, in the limit 
Ri n — » Rout, we have an Einstein thin shell of counter rotating particles [3(, [4|. 

Regardless of the limiting process, we may consider thin shells of counter rotating parti- 
cles and, as is clear after the appropriate equations are set, still be able to impose further 
restrictions on the particle motions. These systems, that are among the simplest dynamical 
models of matter in spherically symmetric space-times, were studied thirty years ago by 
Evans [5j, who considered a manifold with an embedded time like hypersurface, also with 
spherical symmetry, that represents the world tube of a thin shell. All matter content is 
within the hypersurface so that it is singular and the Darmois-Israel formalism [6( is appli- 
cable. In order to get definite equations of motion he also made two important assumptions 
that we will discuss later: i. particles inside the shell move along geodesies of the induced 
metric, and ii. the absolute value of the angular momentum is the same for all particles. 

Models like that of Evans (as special cases of application of singular hypersurface theory) 
have been studied by a number of authors in different applications, including, for example, 
the analysis of certain properties of the evolution of star clusters [7|. The stability of single 
shells under some general assumptions on the type of matter contents has also been consid- 
ered in the literature [8( but, as far as we know, the stability of the trajectories followed by 
particles within the shell, or the stability of a shell under possible splitting into separate com- 
ponents have not been analyzed. This is important, because, as we shall show, the equations 
of motion are satisfied also for shells that are intrinsically unstable, and therefore, certain 
results regarding the maximum compactness (ratio of radius to mass) may need revision 
when stability is taken in consideration. In the case of shells undergoing dynamical evolu- 
tion, a new feature related to stability that results from our analysis is the possibility that, 
under evolution, an initially stable single shell may split into several separate components. 

The plan of this paper is as follows. In the next Section we review the formulation of 
the model. In Section III we obtain the equations of motion for a shell with an arbitrary 
distribution of angular momentum for the constituent particles. Then, in Section IV, we 
consider some examples of static shells, including a two component shell, that shows explic- 
itly that the equations of motion do not provide un upper bound to the possible angular 
momentum of the particles. In Section V we analyze single particle orbits, and obtain gen- 
eral conditions for the stability of the shell under single particle "evaporation" . In Section 
VI we review the stability of the shells as a whole, finding results in agreement with those 
of previous authors. Sections VII and VIII are devoted to the analysis of the stability under 
shell division of both static and dynamic two component shells. Explicit stability conditions 
that relate the angular momenta to the radius and masses are obtained, and some explicit 



examples are analyzed. Although restricted to a two component shell, the results should be 
qualitatively applicable to more general configurations. We also include an Appendix where 
we show that the thin shell limit of an Einstein shell is a shell with a unique continuous 
distribution of angular momenta. We close with some additional comments on the results 
obtained. 



II. THIN SHELLS DYNAMICS IN GENERAL RELATIVITY 

To construct a general model of a thin, spherically symmetric shell of matter we consider 
a spherically symmetric space-time in which the thin shell corresponds to a timelike 3-surface 
£ embedded in this space-time. S is the common boundary of two disjoint manifolds: Aii 
and M.11. Spherical symmetry implies that the line elements on Aii and A4n can be written 
as: 

dsj n = -F^udt] n + — — drj u + rj n d^ 2 (1) 
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where F n = 1 - 2M H /r u and Fj = 1 - 2M I /r I , and dtt = d9 2 + sin 2 
We write the metric induced on X as 

dsl = dr 2 + R{r) 2 dQ (2) 

where (r,9,<f)) are coordinates on the shell, with r the proper time of an observer at fixed 
(9, (ft) on the shell. We take 77 < R(t), and r n > R{t). Therefore S can be described by: 
Si,n(rijj, tiji, 9, 0) = rjjj — R(t) = 0, and R(r) is the radius of the shell. 

We will apply the Darmois-Israel formalism for singular hyper surfaces to obtain the 
equations of motion of the shell. The non vanishing components of the extrinsic curvatures 
of S, respectively in M.i and Ain, in (t,9,4>) coordinates are: 



A". 



lu R 2 R + M IyII 



TT 



R^R(RR 2 + R-2M ItII ) 
K\f = -tJr(RR 2 + R-2M iji ) (3) 
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where R = dR/dr, and R = d 2 Rjdr 2 . 

We can write the energy-momentum tensor for a singular hypersurface S as: 

T^(x a ) = 5(S(x a ))S^ (4) 

Then, Israel's junction conditions |6( are given by: 

(n a n a )([K l:i } - 3 gij [tr(K)]) = nS^efe" = kS %] (5) 

where for a tensor A, we define [A] = A 11 — A 1 , e; stands for d/dr, d/d9 or d/d(p on the 
shell, and n a is the unit normal vector of the shell. 
Applying equation (Q we have: 
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These equations are completely general, and, once they are supplemented with an equa- 
tion of state for the matter contents, they completely determine the dynamics of the shell. 
A very simple example is given by a shell of "dust", where, Sg = S^ = 0. Then, (JSJ) alone, 
up to integration constants, determines the dynamics. 

In this paper we will be interested in the case where the matter comprising the shell 
is made up of non interacting particles, each particle being subject only to the (average) 
gravitational field generated by the rest, and, possibly by some other external source. In 
principle, then, one would expect that the particles move along geodesies of the total grav- 
itational field. We notice, however, that the Christoffel symbols are discontinuous on the 
3-surface. In particular, a world line contained in the shell that is geodesic with respect to 
the connection of region I is not a geodesic for the connection of region II, and viceversa, 
and more generally, a geodesic relative to the induced metric on the shell does not corre- 
spond to a geodesic in either A4j or M.n. We must, therefore, make more precise what we 
mean by "non interacting particles" in this context, since we expect the world lines of such 
particles to be geodesies in some well defined sense. [31] We will make the same assumption 
as Evans J5(: particles inside the shell move along geodesies of the induced metric. As a 
consequence of this assumption, angular momentum is a constant of the motion for each 
particle. Moreover, the restriction to spherically symmetric configurations, implies that for 
every particle moving with a certain velocity in a certain direction relative to the shell we 
must have another moving in the opposite direction, with the same speed, i.e., "counter - 
rotating" particles, as in the Einstein model [2j, or, for the dynamical case, as in [5J. We 
notice that this assumption does not restrict the possible values of the particle's angular 
momentum to a single value. A particular example of a shell with a continuous distribution 
of values of the particle's angular momenta, is given by the limit of shells of counter rotating 
particles with non vanishing thickness (Einstein shells) as the thickness goes to zero. This is 
discussed in more detail in the Appendix, where we show that the limit thin shell must con- 
tain particles with a range of values of angular momentum, so that, in particular, the model 
analyzed by Evans does not correspond to the limit of an Einstein shell. It is, nevertheless, 
expected that different limiting distributions, either discrete or continuous, would result if 
instead of an Einstein shell one considers more general shells of Vlasov type matter. [9( 

Going back to the general formulation of the problem, and taking into account the sym- 
metries of the problem, the 4-velocity of a particle in the shell is completely determined by 
the modulus of its angular momentum per unit rest mass L and the angle x — tan -1 (uq/u</,), 
|5j, and may be written in the form: 
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where \i — T,n, 6, 4>, with n a Gaussian normal coordinate with respect to the shell. There- 
fore, we can write the surface stress-energy tensor as: 



S ij (r) = ^J a(L, r)«*(x, L)u j ( X , L)d X dL (9) 

where i = r,9, <fr, and cr(L, r) is the proper mass density of particles with modulus of angular 
momentum per unit rest mass L, at time r. As we indicate in the next Section, particle 
number conservation leads to a more explicit expression for the functional form of cr(L, r). 



III. EQUATIONS OF MOTION 

The treatment in [5j was restricted to particles with a single value of L. This considerably 
simplifies the expression for S 1 *- 5 and, as a consequence, the form of the equations of motion 
for the shell. Here we consider the more general case, where the angular momentum of the 
particles can take arbitrary values. 

We can use now particle number conservation to obtain first integrals for the Einstein 
equations. Equation (3.17) of 10j implies that the integral on a 2-surface contained in the 



shell world tube defined by r = t : 

- / a{L,r)<ni a u a dL (10) 

A 4 Jt=tq 

where \x is the proper mass of the individual particles, is independent of the choice of r 
and represents the total number of particles inside the shell. Moreover, because of the lack 
of interaction between the shell particles, the number density of particles whose angular 
momentum is L: 

n(L) = - J a(L, r)d£X = - J a(L, r„) (l + ^yj ) R ( T o? sm(0)d8d((> 

= —a(L,T )R(T WL* + R(T )* (11) 
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is also constant. Taking this into account, the components S^ can be written: 



** = i^m dkg 



A simple consequence of I u2l) is that 
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^RJ VL 2 + R 2 ~ V ; 

and, therefore, 

- S T T > 2S e e (14) 

with equality approached only in the limit L 2 — > oo, corresponding to ultra relativistic (or 
massless) particles. If we define a = Sq 9 /(—S T T ), then, 

n{L)L2 dL, (15) 



Jn(L)VR 2 + L 2 dLj 2VL 2 + R 2 



and equation (TT4|) implies < a < 1/2, and from ([7]), we have 
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Notice that, given n(L), a depends only on R, and, therefore, ([TBI) is an integro differential 
equation for R(r). 

Equation (fT2"j) together with (J6J can be used to obtain: 
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and it can be checked that (ITTjl is a first integral of (1161) . 

Defining /(.R) = / n(L)\/R 2 + L 2 dL, this expression can be written as (C = 2/^): 
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Equation (fTHj) is a generalization of equation (2.35) in [5J for an arbitrary angular mo- 
mentum distribution n(L). We can see that the qualitative features of the motion of the 
shell are the same as those described by Evans. V(R), f(R) and its first derivative have the 
following asymptotic behavior: 
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where N is the total number of particles and (L) is the mean value of L. The "binding 
energy per particle" is the constant -1 + A(M n - M 7 ) 2 /(C 2 iV 2 ). If 2(M H - Mi) > CN 
the binding energy is positive and the shell can reach infinity, otherwise, the shell may only 
reach a maximum possible radius. There are also solutions where the shell has a lower 
turning point, that is, the radius reaches a minimum value, which can be obtained solving 
V(R) = 0, for R. Some care should be exercised here, however, because (TTB1 does not take 
automatically into account the signs of the square roots in (fTTI) . In the following Sections 
we consider some simple applications of this formalism, paying particular attention to the 
dynamical stability of the resulting models. 



IV. STATIC SHELLS 

The simplest dynamics results, of course, if the shell is static, with a fixed radius R = Rq. 
In this case we must have both R = 0, and R = 0. From (TT81) . this implies both V\r = r = 0, 



and (dV/dR)\n = RQ = 0, and, therefore, the static shell radius corresponds to a critical point 
of V(R). These conditions are equivalent to setting R = 0, and R = in (jSj) and (J7J), i.e., 



2 
R 

R-M, 



1- 



2M 



lz 



/? 



1- 



2M 7 



~ [n{L)VR 2 + L*dL (21) 



ri 



fl-Mj 



R^R{R - 2M n ) R^/R{R - 2M 7 



2,U 

^3 



n{L)L 2 

VR 2 + L 2 ' 



-AL 



(22) 



For a given distribution function n(L) these equations relate the equilibrium radius R , 
and the masses Mj, and Mjj with the parameters that characterize n(L). Some simple 
examples considered below illustrate this point. 



A. Single Component Shell 

Static single component shells are a particular instance of the systems analyzed in [5|. In 
this case we have n(L) = N5(L — L ), that is, all particles have the same angular momentum 
Lq. The resulting equations are: 
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R-M! 



R 3 

CNLl 



R^R{R - 2M n ) R^R{R - 2Mj) 2i?V i?2 + L o 
It is interesting to solve these equations for L\ and CN. We find, 
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3VR - 2M IV /R - 2M U - R 

For reasons that will be made clear below, we shall call this expression for Lq, the critical 
value of the angular momentum for a shell of radius R, with exterior mass Mu, and interior 
mass Mj. For CN we get, 
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(y/R - 2Mi - y/R - 2M II )Ry/Qy/R~^2M~ I y/R - 2M7i - 2R 



(26) 



^R - 2M IV /R - 2M U 

and we must require Mu > Mj > 0, so that CN > 0. We see that meaningful values of Lq 
and CN are only obtained provided 3^-R — 2Mjy/R — 2Mjj > R, and this implies 
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where (3 = Mj/Mjj and the lower bound on the right hand side of (J2"T|) corresponds to j3 
In this case (Mj = 0), we also have, 

R 2 






3M H -R+ y/R(R - 2M n ) 
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The bound given by (1271) has been considered as the correct lower bound for the com- 
pactness (i.e. ratio R/Mu) for this type of thin shells, with implications, for instance, for 
the maximum gravitational red shift for light emitted from their surface |9j . However, as we 
shall prove in the following sections in our analysis of the stability of the shells, this bound 
cannot be approached if the shell dynamics is governed only by gravitational interactions. 



B. Shells with two components 

A detailed analysis of the dynamics of shells containing particles with more than one 
value of L is rather complex in general. Nevertheless, even if we restrict to two components, 
that is, to two values of L, some interesting features appear that deserve closer examination. 
In this case we have that the number density for particles with angular momentum L can 
be written as n(L) = Ni5(L — L{) + N 2 5(L — L 2 ), and, therefore, after some simplifications, 
we have, 



2M f r 2Mn CN 1 ^R 2 + L 2 CN 2 ^R 2 + L 2 
R V R 2R 2 + 2R 2 



(29) 
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We may solve these equations for Aq, and N 2 , 
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where Lq is the critical value of L, given by (1251) . This result implies that Aq and A*2 will be 
positive, as required by the particle interpretation, provided only that either L 2 > L > L\, 
or L\ > L > L 2 , that is, that one the Lj's is larger and the other smaller than the 
critical value L , for the given values of R and M. In fact, if we consider a static shell 
with any number of components Lj, it is not difficult to show that if the correspond- 
ing Ni are positive, then there must be values both larger and smaller than L$ among the Lj. 

This can be seen as follows. For a static shell with n components, we have Ni, i — 1 . . . n, 
and corresponding Li, i — 1 . . . n, and we may assume L\ < L 2 < ■ ■ ■ < L n . Then from the 
equations generalizing (129]) . and (130~1) . and the definition of L , we have, 



CN^R 2 + L 2 = J2 CN i\/ R2 + L l ( 32 ) 

T 2 n r 2 

CN - ° = YCN^^L= (33) 



Combining (I3"2"j) . and (1321) . we have, 



R 2 n F? 2 

CN - = V CN , (34) 



then, since all A/j > 0, if we assume that L\ > Lq, from (l32l we have 



CNJR 2 + L 2 > Jr 2 + L 2 J2 CN (35) 
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while from (1341) we obtain 
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which contradicts (1351) . and, therefore, the assumption that L\ > L , and Aj > 0. Similarly, 
we may prove that we must have L n > Lo, with the equal sign holding only if Lj = Lq for 
all 2. 

Going back now to the two component shell, we find that (}3"Tj) is a very intriguing result, 
because it implies that the angular momentum of the particles in one of the components, 
say L 2 , can take arbitrarily large values for any admissible R and M, and we still have an 
admissible solution. But this is also true even if R » M, where the regime is essentially 
Newtonian, and since large L 2 implies large speeds for the particles, it is hard to see how 
the particles would remain bound to the shell, if only the gravitational attraction of the 
other particles is responsible for the particle trajectory. Similar problems arise when we 
consider more than two components. As a first step towards understanding this property 
of the solutions of the equations of motion of the shell, in the next Section we analyze the 
stability of single particle orbits. 

V. STABILITY OF SINGLE PARTICLE ORBITS 

We recall that the hypersurface S containing the shell may be considered as a boundary 
for either Aij or M.u. Let us fix now our attention on one of the particles moving on the 
shell with 4- velocity u a . From metric continuity, we may consider this particle as contained 
in either A4i or A4n, with the same 4- velocity u a . In general, from the point of view of A4i 
or M.n-1 namely considering its world line as a curve on either part of the space time, the 
particle does not follow a geodesic path, but rather, it has a non vanishing 4-acceleration. 
To analyze the stability of the orbit of a single particle moving on the shell with 4- velocity 
u a , we consider an infinitesimal (radial) displacement of its world line, so that the particle 
may now be considered as moving freely in either Aii or Aiu, depending on the direction of 
the displacement. Because it moves now along a geodesic close, but outside E, the particle 
will acquire an acceleration relative to the shell. The motion will therefore be stable if the 
acceleration points towards the shell, but it will be unstable is the relative acceleration points 
away from the shell. To illustrate this point, and before we discuss the general case, let us 
consider a static shell, of radius R and (external) mass Mu, which, in accordance with our 
previous discussion may contain particles with arbitrarily large angular momentum L p . The 
metric in the external region A4u, (r > R) may be written as, 

ds 2 = -(1 - 2M H /r)dt 2 + (1 - 2M II /r)- 1 dr 2 + dr 2 dQ 2 (37) 

Assuming equatorial motion, (9 = 7r/2), and calling a the proper time along the particle 
world line, the geodesic equations imply dtp/ da = L p /r 2 , 

d±Y r 2 (dr\ 2 L 2 p + r 2 

da J (r - 2M n ) 2 \da J r(r-2M n ) K ' 



and, 



d 2 r M n r-3M n , . , 

ss. = -7? + ^r^i? (39) 
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Since, for a particle on the shell (r = R) we have drjda = 0, for R > 3Mn the motion 
will be stable against "single particle evaporation", only if, for R > 3Mjj, we have L 2 P < 
M n R 2 /(R—3M II ), while for R < 3M n , all orbits are stable. This is, of course, a consequence 
of the fact, that in Schwarzschild's space time, time like geodesies cannot have a lower turning 
for r < 3M. We also remark that for a single component shell, L p = Lq, and, as can be 
checked, the shell is always stable under single particle evaporation. 

To analyze the more general case of a dynamic shell it will be advantageous to introduce 
a coordinate system (r, p, 9, </>), adapted to the shell, where p is a gaussian coordinate normal 
to S. We are only interested in the limit p —* 0. Then, keeping only terms up to linear 
order in p, in the neighbourhood of the shell (p = 0), we may write the metric in the form, 

ds 2 = -(1 + Ap)dr 2 + 2Bpdrdp + (1 + Cp)dp 2 + (R 2 + Dp)dn 2 + 0(p 2 ) (40) 

where A, B, C and D are functions of r given, 

2(R 2 R + M±) 



B 



R^J(RR 2 + R-2M±)R 
(6Af± - 3R-2RR 2 ) M±R 



R(R - 2M± 



1 2 



2M±(2M± -R- 2RR 2 )\ R(RR 2 + R- 2M±) 

R 2 {R-2M ± ) 2 y ' 



D = 2^/(RR 2 + R-2M ± )R 

where M± = Mjj for p > (the exterior region), and M± = M/ for p < (the interior 
region). R = R(t) is the shell radius, and dots indicate derivatives with respect to r. We 
also require Mu > Mi, for consistency. This metric is continuous, with discontinuous first 
derivatives. For a particle instantaneously moving on the shell we have dp/dt] = 0, where r] 
is the proper time along the particle world line. Considering again an infinitesimal (positive) 
radial displacement, so that the particle follows a geodesic path in A4n, the motion will be 
stable only if for the resulting world line d 2 pjdif < 0. Without loss of generality we may 
consider equatorial motion with 9 = ir/2. Then we have, 

u^ = ^ 
R 2 

(u T ) 2 = 1 + ^ (42) 

where u a = dx a jdr\. Then, the geodesic equations of motion imply, 

d 2 p Mu + R 2 R 2 R- 3M H - R 2 R + RR 2 

d7]2 rJr(RR 2 + R- 2Mn) R^R(RR 2 + R-2M n ) 

This results reduces to fl39l) for the static case R = 0, R = 0, the apparent difference being 
due to the different definitions of r and p. For R > 3Mn and R < 0, this result implies that 
for sufficiently large L p , we have d 2 p/dr] 2 > 0, and the resulting motion is unstable. The, 
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perhaps, unexpected result is that for sufficiently large R > 0, that is, a shell undergoing 
fast accelerated expansion, we may have d 2 p/dr] 2 < for all L p , so that the shell is stable 
against "single particle evaporation". One can check that this condition is consistent with 
the assumption that the shell is made out of counter rotating particles, but we shall not give 
the details here. 

We may also consider a negative infinitesimal displacement. In this case, we obtain again 
(|43p . but now with Mj, replacing Mjj, and now stability corresponds to dPp/drj 1 < 0. The 
analysis is similar to that for (T431) . and will not be repeated here. 

The foregoing considerations apply only to what we might call " single particle evap- 
oration". Perhaps of much greater interest is the stability of the shell as a whole. Here, 
however, we must distinguish different modes. We shall restrict to modes that preserve the 
spherical symmetry, and consider two cases. The first, analyzed in the next Section, will 
be the stability of stationary configurations, and the second, in the case of shells with a 
distribution of values of L, the stability under separation of components. 

VI. STABILITY OF STATIC SOLUTIONS 

To analyze the stability of static solutions we consider again ([TBI) . Setting now R = R = 0, 
we find, 

1 







2^o 



IMi / 1M U 



(44) 



where Rn is the static shell radius. This can be solved for a, 



.„. R - VRo - 2Mj^R - 2M n ,,_. 

a(R Q ) = . - . = — (45) 

Ay/Ro - 2M iy /Ro - 2M/J V ; 



Since a(R ) < 1/2, we obtain the bound [9| 



9,_ ., , .3 rrr, — ttt^~ ttt: . 9 



Ro > g (M/ + M n ) + - V 9 (M 7 - M n y + 4MjMjj > -M n (46) 

for any shell of counter rotating particles. The lowest value can only be approached for 
Mj = 0. This is the same as ( 1271) . but in this case there is no restriction on n(L). 

We consider now a small perturbation of the shell radius that preserves spherical symme- 
try and the angular momentum distribution n(L). Then a may be considered as a function 
of R only, and we have a' = da/dR < 0. We introduce 5R(t) such that R = Rq + 5R(t), and 
expand (Tl6l) to first order in SR and its time derivatives. After some simplifications using 
also (145]) . we find, 

d 2 5R 2VR - 2M /v /#o - 2Mjj „„..„_, (Mj + M n )R - AMjM n 

^ = Rl a {Ro)SR + 2RI(R -2M I )(R -2M II ) SR (4?) 

Since we take Mjj > Mj, and we must have Rq > 2Mji, the first term on the right hand 
side of (1471) is always negative, contributing to stability, while the second term is always 
positive, tending to make the system unstable. 
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We may consider now two limits. First, from (fl~5l) . for large i?o, assuming R » L in 
T5l). we have, 



fn(L)L 2 dL 1 2 

o' ~ — — - — = < L > 

" 2R 2 fn(L)dL 2R 2 

f n(L)L 2 dL 1 o . , 

a <*»> - -ifjiiMi = "^ < L " (48) 

where angle brackets indicate average with respect to n(L). Then, from the the staticity 
condition (1451) . we have, 

< L 2 >~ - (Mi + Mn) Rq (49) 



and, replacing in (J4"T|) . we find, 

d 2 5i? (M 7 + Mn) 



5R (50) 



d£ 2 2i?g 

Therefore, all static shells are stable (as regards the mode considered here) for sufficiently 
large i? - 

The other limit is as Rq approaches the minimum value allowed by (j4"5l) and a < 1/2. 
But this requires a ~ 1/2, which in turn requires large values of L as compared with i? in 
(TT5|) . In this case we have, 



1 R 2 

a - 



f n{L)jdL 



2 2fn(L)LdL 



1 


R 2 


<L~ 


x > 


2 
< 


2 
L- 1 


<L 

> 


> 



■ R °^TT (51) 



Then, as i?o approaches the minimum value, we have a'(Ro) — ► 0, and, therefore, from (j4"Tj) . 
we conclude that a// shells are unstable in this limit. To find the minimum stable radius R s , 
from (J4"Tj) . we see that for stability we must have, 

7/n< (M 7 + M n ) Rs ~ ±MiMn f . 

{ s) ~ A(R S - 2Mif/ 2 (R s - 2M H y/ 2 { ' 

from which we may obtain the minimum equilibrium radius R s , for a given distribution 
function n(L). As an example, for a single component shell we have a(R s ) = Lq/(2(Lq+R 2 )), 
and f|45|) implies, 

l2 _ R 2 S {Rs - VRs - 2M^R S - 2M n ) 
3y/R s - 2M iy /R 8 - 2M U - R s ^ ' 



Then, from (1521) we find that R s is a solution of, 

AR 2 - 9 (Mi + Mn) R s + 20M J7 M 7 2i? s 2 - 3 (M 7 + M /7 ) i? s + 6M /7 M 7 
2y/R a - 2M^R S - 1M n Rs 



(54) 



The appropriate root of (1541) grows smoothly as Mi grows from Mi = up to Mi = Mn. 
Near those extremes it satisfies, 



. 51 3V33 \ / 3 19V33 . 

' 16 + ^6~ ) MlI+ 16 + ^8"I Ul ' 4 - 26M »- 
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for Mj — > 0, and, 

R a ~ 6M„ - 3 (M n - Mi) (56) 

for Mi — > Mjj. Therefore, the lower bound on R/Mjj set by the stability criterion is almost 
twice as large as that given by ( |27|) . (See also [81] for related results.) 



VII. STABILITY UNDER SHELL DIVISION OF STATIC TWO COMPONENT 

SHELLS. 

As indicated in Section IV, in the case of shells with a distribution of values of L, one can 
have solutions of the equations of motion even if the shell contains particles with arbitrarily 
large values of L. As shown in Section V, this may lead to instabilities under "single particle 
evaporation" , but it is natural to question the possible stability of such configurations as 
a whole. It turns out that the general problem is rather cumbersome and hard to handle. 
Nevertheless, one can get important information and insight by considering first static shells. 
In this Section we consider the stability of such shells under the further, but not trivial, 
simplification of assuming only two components, one with angular momentum L\, and other 
with L 2 , with L 2 > L\. This configuration was already analyzed in Section IV, to obtain the 
static solution, assuming an inner mass Mj, outer mass Mjj, and shell radius R$. The central 
idea here will be to assume that the L\ and L 2 components are infinitesimally displaced from 
the static radius Ro, the first to the inside and the second to the outside, so that we have 
now an intermediate vacuum region separating the two resulting shells. Let Mj ni be the 
mass parameter for this intermediate region. On account of the derivations in Section III, 
and assuming that the shells acquire no velocity in the displacement, we have, 



d 2 R 1 


1 
2-Ro 

1 
2-Ro 


1- 

1- 


( 2L? \A 


t- 


2M int 


i- 


2Mj 
R 


(57) 


drl 


\m+Li' ) 
( 2Ll 1 1\ 


Ro \ 


d 2 R 2 


iA- 


2M int 


/- 


2M n 
R 


(58) 


drl 


V«S + ii ) 


Ro \ 



where Ri and R 2 are the shell radii and T\ and t 2 the corresponding proper times. We also 
have, 



/ 


2M int 
Ro 


I 2Mi CN^Rl + Ll 
V ^o 2R 2 


/ 


2M int 


CN 2y /R 2 + Lj I 2M H 
2R 2 V Ro 


R 



(59) 
(60) 



We may use now (I3~T1) to eliminate CN\ and CN 2 from these equations and solve for M int . 
We find, 



1- 



2Mi 



int 



Ro 



{Ll-L\){R 



{LI 



L\){Rl 



LI) 



LI) 



{Ll-Ll){Rl + L\) 



{LI 



L\){Rl 



LI) 



1 



2Mj 

Ro 



2M 
Ro 



ii 



(61) 
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which is of the form, 



2M int 2M T 1M n 

-^ =w i 1 -i^ + H 1 -^- (62) 

with Wi > and w\ + u> 2 = 1 , so that we always have Mj < Mi nt < Mjj, as expected. 

From these equations we immediately conclude that the shell will be stable under separa- 
tion into two single component shells only if d 2 Ri/dr 2 > and d 2 R2Jdr 2 < O5 but unstable 
otherwise. The critical condition is then d 2 Ri/dr 2 = 0. It is interesting that both conditions 
are satisfied if L\ and L 2 are related by, 

(3L 2 + R 2 )VR -2M n (3L 2 + R 2 )VR -2M I 

(R 2 + L 2 ) (R 2 + L{) [M) 

We can check that this is indeed a critical condition by computing, e.g., the partial 
derivative of the R.H.S. of (158]) with respect to L 2 , keeping L 1; i? , Mj, and M n fixed, 
evaluated assuming f[6"3"l) holds. The result is, 

4L 2 # 2 (L 2 2 - L 2 ) (R - 2 M H ) ^R -2M! 

— — ^^ 7) > (64) 

(VRo-2M T - VRo-lMu) (3 L 2 2 + .Ro 2 ) (R 2 + L 2 2 ) (R 2 + L 2 ) 

Therefore, the shell becomes unstable for values of L 2 larger than the one satisfying (1631) . 
Similar results hold for L 1; except that here the shell becomes unstable for values of L\ less 
than the critical one given by (163]) . 

The relation (1631) must also comply with the conditions L\ < L 2 ,, and L\<L\. We can 
see from (16"3"1) that, for L 2 > 0, L 2 and L 2 are monotonic functions of each other. Then, the 
lowest possible value for stability for , e.g., L\ is, 

(L 2 ) . = MlR2 ° (65) 

V l/mm R Q -3M! y ' 

because for this value we have L\ = Lq, the minimum possible value for L 2 . The value 
(|65|) coincides with the single particle stability limit, but we may also check that for this 
configuration, with L\ = Lq, we have CNi = 0, and the particles with angular momentum 
equal to L\ become test particles. Similarly, the highest possible value for L\ is L\ = L . 
Here stability requires, 

(L 2 ) = MnR2 ° (66) 

which is again the single particle stability limit, but for L 2 , with similar considerations as 
for L\. The important general conclusion from this analysis is that the criterion for stability 
under shell division is stronger than that for single particle evaporation, since it leads to 
critical Lj values that are either lower ion L 2 , or higher than L\ than those required from 
single particle stability. 

It is interesting to check also that a static shell with a single component is stable under 
splitting. We consider again Eqs. (I57II58II59||6"U|) . with L\ = L 2 = L . From (I59|I60I) we find, 



2M int N, 2M n N 2 2Mj 

R iV! + iV 2 V R N, + N 2 \ R [ } 
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with Ni, and N 2 restricted by, 
C(N t + N 2 ) -- 



2RI 



vw^n 



2M 1 
R 



1- 



2M 



11 



Rq 



(68) 



but otherwise the ratio N 2 /Ni is arbitrary. Replacing (IBTl) . and L 2 = L in, e.g., fl5E\i 

d 2 fl 2 __ iVi(Vi?o - 2M J7 - ViZp - 2M 7 ) 
drl ~ 2(N ± + N 2 )R Q ^R - 2M 7 

so that the shell is stable for all ratios N\/N 2 . 



<0 



(69) 



VIII. STABILITY UNDER SHELL DIVISION OF NON STATIC TWO 

COMPONENT SHELLS. 



As indicated in the previous Section, the stability analysis in the general non static case 
is, even for a two component shell, considerably more complicated in detail than in the static 
case. The analysis may be carried out along lines similar to those used for the stability under 
single particle evaporation. We first notice that for a two component shell the equations of 
motion are, 



d 2 R 

Ifr 2 



1 

m 



1 + R 2 - (4a 



2M 7 



R 2 \ 1 



2M, 



n 



R 



R 2 



where, 



and, 



a 



2(N 1 ^/R 2 + L 2 + N 2 ^/R 2 + L 



N^ 2 



N 2 L 2 



^R 2 + L\ ^WTL 



(70) 



(71) 



_^Mi + R 2_ A 2M " i p 2 = CN 1 ^R 2 + L 2 t CN 2 ^WTI 



R 



R 



2R 2 



2R 2 



(72) 



where R = R{t), R = dR/dr, and r is the proper time for an observer that moves radially 
with the shell. Notice that (!72j) is actually a first integral of (1701) . A particular solution of 
the equations of motion is determined by fixing, for some particular value of r, say r = To, 
appropriate values of Afj, M n , (M n > Mi), L u L 2 , (L 2 > L x ), CN X > 0, and CN 2 > 0, and 
a corresponding radius R(tq) = Ro, such that ( 1721) may be solved for a real value of dR/dr. 
We consider now, as in the previous Section, for the same value of r = To, an infinitesimal 
displacement of the components, so that we have two shells, one made out of the particles 
with Li, and the other with L 2 , separated by an intermediate empty region with mass M int . 
We assume that Mj, Mjj, L\, L 2 , CNi, are CN 2 , are not modified by this displacement, and 
that the shells have initially vanishing relative velocity. If we look now at the L 2 component, 
its equations of motion should be, 



d 2 R 2 
~d4 



2R 2 



Ri 



(4a 2 



1- 



2M 



hit 



Ri 



Ri 



2M n ■ 2 



(73) 



where, 



and, 



2M; 
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R>. 



Bo. 



a 2 



1 



r2 

^2 



2(R 2 



2Mi 



int. 



R, 



L 2 ) 



+ R 2 + 



CN 2y /B% + L\ 
2R\ 
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(74) 



(75) 



where R 2 = R 2 (r 2 ), and r 2 is now the proper time for the L 2 shell. We assume that for 
r = r we also have r 2 = r , so that R(r = r ) = R 2 {r 2 = To). From the above assumptions 
and metric continuity, we also have dR 2 {r 2 ) / dr 2 = dR(r)/dr. Then, we may use (T75"|) to 
compute M inU and then obtain d 2 R 2 {r 2 )/dr 2 from (I75|) . We might be tempted to conclude 
that the condition for stability of the shell should then be d 2 R 2 (T 2 )/dr 2 < d 2 R(r)/dr 2 , but 
this could in principle be incorrect, because of the different meanings of the proper times r, 
and t 2 . It turns out, however, that this is the correct criterion. We may prove this as follows. 
From the point of view of an observer in the region Mji, the motion of the shell is described 
by a function R = R(tn), and we have the following relations between the derivatives with 
respect to t n and to proper time r (or r 2 ) of the shell, 



dr 



dt 



n 



(R(t) - 2M n ) 



R(t) 2 - 2R(t)M u + R(t) 2 ( 



dR(r) 
dr 



(76) 



from which we may also compute d 2 r/dt 2 - I , and d 2 r 2 jdt 2 u . From these expressions we may 
finally derive, 



d 2 R{t 



ii) 



R 



dt 2 n 



R-2M U \dtn 

M n 



4 ^2; 



dr y d 2 R(r 



R 3 (R - 2M U ) 



dr 2 



R - 2M n + R 



dr 



dt 



ii 



2R - AM n + R 



dr 



dt 



n 



(77) 



and a similar expression with R replaced by R 2 and r by r 2 . But, from the above assump- 
tions, for r = t 2 = To we have dr/dtu = dr 2 /dtn, and then the shell with the larger proper 
radial acceleration will also be seen to have the larger radial acceleration with respect to tn. 
Moreover, the critical condition, i.e., equality of the second derivatives with respect to tn, 
implies also equality of the second derivatives of R, and R 2 , with respect respectively to r 
and t 2 . If we assume that at a certain point, we have R(t) = R{t 2 ), dR/dr = dR 2 /dr 2 , and 
d 2 R/dr 2 = d 2 R 2 /dT 2 , an replace in (1701) to (175|) . we obtain the following condition, 



(3L| + R 2 " 

(R 2 + L 2 ) 






2M 



n 
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dR 

~drr 



(3L 2 + R 
(R 2 + L 2 




1- 



2M T 



R 



(dR 



(78) 



which is seen to be a simple generalization of the critical condition for stability for static 
shells given by fl63l) . It is again remarkable that this is also the critical condition if we 
consider stability for the shell with L\. Since the values of R and dR/dr change as a 
function of r, a shell that is initially stable may become unstable as it expands or contracts, 
and spontaneously split into two separate components, one with angular momentum L\ and 
the other with L 2 . 
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However, for an expanding shell, as R increases, assuming given value of L\, Mj, and 
Ma, dZHD is satisfied for increasing values of L 2 . This can be seen solving (1751) for L\, and 
assuming that R is much larger than the other quantities in (1781) . (we recall that R = dRjdr 
is always bounded, even for R — > oo). The result is, 

L\ = M »- M >R + L\ + M„(M„-M,) + (^) (79) 

2(1 + i? 2 ) X (1 + i? 2 ) 2 V ; V 7 

and, therefore, for sufficiently large R, any L 2 will satisfy the stability condition. Since 
instability results for either L 2 larger (or L\ smaller) than the critical value, we conclude 
that any two component shell necessarily becomes stable if it expands to sufficiently large 
values of R. 

The derivations for the two component shells strongly suggest that these results hold also 
for any multi component shell. If this is correct, a thin shell with a continuous distribution 
of L values might, upon evolution, turn into a thick shell with a continuous plus a singular 



distribution of matter, an interesting result in view of some theorems [11| on the emergence 
of singularities in Vlasov - Einstein systems. We shall consider this point in more detail in 
a separate paper. 

Going back to the stability problem, in principle, given Mi, Mu, L\, L 2 , CN±, and 
CN 2 , we may solve ( J72l) for (dR/dr) 2 , and replace in (178]) to obtain an equation for R in 
terms of these constants, whose solutions contain the values of R where the shell becomes 
unstable. Unfortunately, the explicit expressions are too complicated to allow for a simple 
interpretation of their meaning. Instead, we consider several explicit examples where we 
integrate (1701) . and use the resulting R(t), and dR/r to compute M int , and then the right 
hand side of (1771) . also as a function of r. The shell will become unstable if this is larger 
than d 2 R/dr 2 in (1701) . Similarly with regards to the shell with L\, except that here the 
shell becomes unstable when the corresponding acceleration is smaller than d 2 R/dr 2 , but 
we must emphasize that, as shown above, the instability condition is always simultaneously 
satisfied for both component shells. 

As a first example we consider a solution of the shell equations that corresponds to 
a bounded periodic motion. For this example we chose Mi = 0.45, Mu — 1, C — 1, 
U. = 2.8,L 2 = 3.7, JVi = 0.3119, and N 2 = 0.8044. Figure 1 shows R(r) as a function of r 
for this choice of parameters. We have two turning points, one for R = 8., and the other 
for R = 40.0.... The period is At ~ 940. We have chosen r = at the lower turning point 
R = 8. 

To analyze the stability under shell splitting, we computed d 2 R 2 /dr 2 , and d 2 Ri/drf, as 
functions of r as indicated above, and compared their values with those of d 2 R/dr 2 . The 
relevant results are indicated in Figure 2, where curve a corresponds to d 2 R/dr 2 , curve b to 
d 2 R 2 /dr 2 , and c to d 2 R\/dr 2 . The curve is restricted to 50 < r < 140 to show the point 
where the shell turns from a stable to an unstable motion, at the critical value r ~ 86, 
corresponding to R ~ 17.. The motion is unstable to the left and stable to the right of 
this point. Given the fact that the motion of the shell as a whole is periodic and invariant 
under r — > — r, we conclude that the motion is stable under shell splitting for R >~ 17 and 
unstable for R <~ 17. We must remark, however, that the if the shell starts its motion in 
the stable region, and a splitting does occur at the critical point, the ensuing motion will no 
longer be periodic, because the resulting separate shells evolve with different proper times, 
and, in general, they will have a non vanishing relative velocity when they cross each other 



again. This kind of motions has been investigated by Barkov, et.al., [7] who find a very 
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FIG. 1: R(t) as a function of r for Mi = 0.45, M H = 1, 
and N2 = 0.8044. The motion is periodic and bounded. 



C = 1, Li = 2.8,L 2 = 3.7, JVi 



0.3119, 



complex chaotic like pattern for the motions of the separate shells. Thus, in this example, 
a highly ordered initial state, corresponding to a single shell, may turn into a very complex, 
non periodic motion of the separate shells. The analysis of the motion of the shells after 
separation is, however, outside the scope of the present analysis. 

In general, it is also possible to find choices of parameters such that, for bounded periodic 
motions, the shell is unstable at all points, or, on the contrary, it is stable at all points. 

As a second example we take M 7 = 0.6295, M u = 1, C = 1, L ± = 10, L 2 = 11, Ai = 0.2, 
and 7V2 = 0.3. Figure 3 shows R{t) as a function of r for this choice of parameters. We have 
a single turning point for R = 8., and we have chosen r = at this point. Here we have two 
critical points, one at r ~ 6.03... corresponding to R(t) ~ 9.10..., (see Fig. 4) and the other 
for r ~ 232.5..., corresponding to R(t) ~ 246.7... (see Fig. 5). The motion is unstable for 
R in this range, but stable near the turning point, and for R larger than the larger critical 
value, in accordance with our previous discussion on the behaviour for large R. 

IX. COMMENTS AND CONCLUSIONS. 



One question that arises in trying to interpret the instabilities we have found in this 
paper is why do we have this ambiguity at the critical points, and unstable regions, where, 
both a single shell and a split pair are compatible with the equations of motion for the 
particles, that are explicitly assumed to be collisionless, and interact only gravitationally. A 
possible answer is that the equations of motion for the shell contain implicitly a "constraint" 
that forces the particles to remain on the shell. The effect of the constraint is irrelevant 
in the stable region, where, at least in the Newtonian limit, one can see a "wedge" shaped 
potential, keeping the particles at the bottom of the wedge, but is crucial in the unstable 
region, where the wedge can no longer restrain the particles. This "constraint" , would have 
to correspond to a massless but both stiff and elastic structure, and , therefore, we consider 
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FIG. 2: Curve a (solid line) corresponds to the acceleration d 2 R/dr 2 of the shell, curve b (dotted 
line), to d 2 R2/dr2, computed as the right hand side of (|73p . and curve c (dashed line), to d 2 Ri/dr 2 , 
similarly computed. Only the region where the shell becomes critical is depicted. The critical value 
of r corresponds to the crossing of the curves, where the three accelerations are equal. To the right 
of this point we have d 2 R2JdT 2 < d 2 R/dr 2 , and d 2 Ri/drf > d 2 R/dr 2 , and the shell is stable. The 
opposite situation holds for r less than the critical value. A vertical unit in this graph corresponds 
to 10" 5 . 
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FIG. 3: R(t) as a function of r for Mj = 0.6295, M H = 
N2 = 0.3. The motion has a lower turning point at R 
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FIG. 4: Curve a (solid line) corresponds to the acceleration d 2 R/dr 2 of the shell, curve b (dotted 
line), to d 2 R2/dr2, computed as the right hand side of (|73p . and curve c (dashed line), to d 2 R\/dr 2 , 
similarly computed, and corresponding to the parameters in Fig. 3. Only the region near the critical 
point at r ~ 6.03... is depicted. The critical value of r corresponds to the crossing of the curves, 
where the three accelerations are equal. To the left of this point we have d 2 R2Jdr 2 < d 2 R/dr 2 , 
and d 2 R\/dr 2 > d 2 R/dr 2 , and the shell is stable. The opposite situation holds for r to the right 
of this critical value (but less than the higher critical value displayed in Fig. 5). A vertical unit in 
this graph corresponds to 10~ 3 . 

the splitting interpretation as the more physical one in our case. This would clearly have to 
be taken in consideration in any application to physically meaningful systems. 

APPENDIX: THE THICK TO THIN SHELL LIMIT FOR STATIC EINSTEIN 

SHELLS 



We consider a static thick Einstein shell, that is, a spherically symmetric space time 
where the matter contents is made out of equal mass, non interacting particles moving 
along circular geodesies, and confined to a spherical shell of non vanishing thickness. The 
metric for this space time may be written in the form, 



ds 2 = -e^dt 2 + 1 



2m(r] 



dr 2 + r 2 (d6 2 + sin 2 6dcj) 2 ) 



(A.l) 



We assume a central source of mass Mi, and that the shell is restricted to Ri < r < R , 
with 2Mi < Ri. Then, for 2Mj < r < R t we have, 



m(r) = Mi ; e v(r) = 1 - 



2Mi 



(A.2) 



Because of the positivity of the energy density, (see, Eq. (IA.9J) ). the function m{r) is 
increasing as a function of r in the interval Ri < r < R , (inside the shell), attaining its 
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FIG. 5: Curve a (dotted line) corresponds to d 2 ^/^! — d 2 R/dr 2 , and curve b (dashed line), to 
d 2 Ri/dT2 — d 2 R/dr 2 , computed as in Fig. 4, and corresponding to the parameters in Fig. 3. Only 
the region near the critical point at r ~ 232.5... is depicted. The critical value of r corresponds 
to the crossing of the curves, where the three accelerations are equal. To the right of this point 
we have d 2 R 2 /dT$ < d 2 R/dr 2 , and d 2 R\ldr\ > d 2 R/dr 2 , and the shell is stable. The opposite 
situation holds for r less than this critical value (but larger than the lower critical value displayed 
in Fig. 4). A vertical unit in this graph corresponds to 10 -7 . 



maximum value m(r) = Mjj, at the outer boundary r = R of the shell. Then, for r > R 
we have, 

2Mii 



m[r\ 



M n ; e"M = C (l - ^l\ 



(A.3) 



where C is a constant. 

The stress-energy-momentum tensor is given by, 

T ab (x a ) = /i _1 n(r)(p a p{,) 



avg 



(A.4) 



where \i is the mass of the particles, n(r) is proportional to the proper particle number 
density, and the particle 4-momenta p a are averaged over all space directions, compatible 
with the condition of circular geodesic motion. Since p r = for all particles, this implies 
that all components of the form T ra vanish. In particular, imposing T r r = on Einstein's 
equations we find, 

dv 2m(r) 



dr r(r — 2m(r)) 
Moreover, the condition of circular orbits implies that for all particles p a satisfies, 



PaP 



fl 



2 t- 2 

_Mr) ( t\ 2 , Z 1 L 



if)' 



(A.5) 



(A.6) 
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where L is the particle's angular momentum per unit mass. A simple computation then 
shows that the only non vanishing components of T^ are, 



T t = -fxn{r) 



1 - — 



-P (A.7) 



Te° = T/ = pj^ (A.8) 

where p is the energy density. From Einstein's equations we have, 

2 dm , . „. 

r A dr 

We shall be interested in the thin shell limit, where Ri — > R — > Rq, where Rq is some 
radius. Then, unless the shell is empty, p approaches a Dirac's 5 form, and so does n(r), 
which is therefore singular in this limit. Nevertheless, as we shall show, if instead of con- 
sidering the number of particles at a given r, we look at the distribution of values of L 2 , we 
find that as the shell becomes thinner and thinner, this distribution approaches a unique 
smooth form, that depends only on Mi, Mn, and, and the radius Rq of the limiting thin 
shell. This can be seen as follows. 

The condition that the world lines of the particles be geodesies of the metric (1A.1I) imposes 
that, 

i 2 = ^^T (A.10) 

r — 6m[r) 

which makes explicit the condition that no part of the shell may have r < 3m(r). In 
particular, we must have R > 3Mjj, so that thick shells cannot be made more compact 
than thin shells. From (1A.10I) we get, 

dL 2 r 3 dm rm(r — 6m) 

dr (r — 3m) 2 dr (r — 3m) 2 

and therefore dL 2 /dr > for r > 6m(r), or dm/dr sufficiently large. We will be interested 
in the case where Ri —>■ R , where dm/dr becomes arbitrarily large, since m(r) grows from 
Mi to Mn in an arbitrarily small interval of r. In this case, we may assume dL 2 /dr > 0, 
and therefore, we have that L 2 grows monotonically in the interval, 

RIM, £L , £ PJM,, 



Ri - 3Mi ~ ~ Ro- 3M n 
Let now n(r)dr be the total number of particles between r and r + dr. Then we have 



,-^r'^j^ 



ATrpr 2 1 ^ =fin(r)\l + — (A.13) 



Since we assume dL /dr > 0, we may introduce, 



n{L 2 )=n{r)^- 2 (A.14) 
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where n(L 2 )dL 2 is the total number of particles in the shell with (squared) angular momen- 
tum in (L 2 , L 2 + dL 2 ). We may solve (1A.10I) for m{r), 

m(r) = ^-^ (A.15) 

and derive with respect to L 2 to get, 

dm r 3 _ / dr \ , . „ nS 

and we remark that dr/dL 2 — > as i?j — > i? . We may now solve ( 1A.13I) for n, and using 
the previous results, we find, 

^ = ^L 2 + R 2 J/ 2 {L 2 + R 2 ) (A - 1?) 

where we have taken the thin shell limit Ri — > R — > Rq, and L 2 is restricted to the range, 

RIM, ^ nlM„ 



R - 3Afj ~ _ R - 3M n 

Thus we see that in the thin shell limit, an Einstein shell contains particles with a unique 
continuous distribution of angular momentum given by h(L 2 ), with L 2 varying continuously 
in the range ( 1A.18I) . 
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